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Abstract 

Generalized parton distribution (GPD) contains rich information of partons in a 
hadron, including transverse profile, and is also non-perturbative information necessary 
in describing a variety of hard processes, such as meson leptoproduction and double 
deeply virtual Compton scattering (DDVCS). In order to unveil non-perturbative as- 
pects of GPD, we study DDVCS at small x in gravitational dual description. Using the 
complex spin j-plane representation of DDVCS amplitude, we show that GPD is well- 
defined and can be extracted from the amplitude even in the strong coupling regime. 
It also turns out that the saddle point value in the j'-plane representation plays an 
important role; there are two phases in the imaginary part of amplitude of DDVCS and 
GPD, depending on relative position of the saddle point and the leading pole in the 
j-plane, and crossover between them is induced by the change of the kinematical vari- 
ables. The saddle point value also directly controls kinematical variable dependence of 
many observables in one of the two phases, and indeed the dependence is qualitatively 
in nice agreement with HERA measurements. Such observation that the gravity dual 
shares basic properties of the real world QCD suggests that information from BFKL 
theory might be used to reduce error in the gravity dual predictions of the form factor 
and of GPD. This article also serves as a brief summery of a preprint larXiv: 1105. 2999 1 



1 Introduction 



AdS/CFT correspondence and its extension to non-conformal theories have been exploited 
for study of non-perturbative aspects of strongly coupled gauge theories. Hadron spectra, 
coupling constants among them and chiral symmetry breaking have been studied intensively 
in the literature by using gravitational dual descriptions with smooth infra-red non-conformal 
geometries. The gravitational dual approach can be used, however, to study not just static 
properties of strongly coupled gauge theories, but also scattering of hadrons. Indeed, string 
theory or dual resonance model was originally constructed to describe scattering of hadrons. 
Qualitative aspects of hadron scattering can be obtained in gravitational dual descriptions, 
if the background geometry (target space) of string theory is chosen properly [H El E]- 

In this article, we will study 2-body to 2-body scattering of a hadron and a virtual photon 
at high energy in gravitational dual descriptions. This process is called double deeply virtual 
Compton scattering (DDVCS). When the final sate photon is on-shell, it is called deeply 
virtual Compton scattering (DVCS), and is accessible in experiments [1]. Because of QCD 
factorization theorem [5], the DVCS or DDVCS amplitude is obtained as a convolution of 
generalized parton distribution (GPD) |6] and a hard kernel, the latter of which can be 
calculated in perturbative QCD. GPD itself (at a certain factorization scale), however, is 
a non-perturbative object in nature, and cannot be calculated in perturbative QCD. Even 
in determining it by using experimental data, its profile needs to be parametrizec l] based 
on proper understanding on non-perturbative dynamics behind confinement. We thus use 
gravitational dual descriptions to extract theoretical understanding on the GPD profile. 

It is not that we just use a well-developed technique to calculate a specific scattering 
amplitude (or GPD) in this article, however. This article clarifies structure of Pomeron "ex- 
change" amplitudes, how to organize them, as well as their field-theory interpretation. We 
find that a saddle point value of the scattering amplitude in complex spin j-plane represen- 
tation is a key concept in organizing Pomeron amplitudes and in understanding kinematical 
variable dependence of the scattering amplitude. Based on this understanding, sharp cross- 
over behavior is expected in the photon-hadron 2-to-2 scattering amplitude in small x limit. 

This article is meant to be a brief summary of reference [8]. To keep this letter short 
enough, we extracted material mainly from §5 of [8] , and only minimum from other sections, 
imagining people in perturbative QCD community as primary readers of this letter. More 
theoretical aspects of the scattering amplitude in gravity dual, as well as more detailed 
account of the materials in this letter, are found in [S]. 



^See [7] for review articles, which also have extensive list of literatures. 
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2 Amplitude in Gravity Dual 



In order to calculate hadr on-virtual photon scattering amplitude in gravitational dual, one 
needs to adopt a certain holographic model. Since the real world QCD turns from weak 
coupling at high energy into strong coupling at infrared, it is desirable to have a holographic 
model that is faithful to string theory where AdS curvature becomes larger than string scale 
toward UV boundary. Such a model becomes even more realistic, if spontaneous chiral 
symmetry breaking is implemented in it. Our primary goal in this article, however, is not in 
pursuing precision in numerical calculation (as lattice QCD does) by setting up a perfectly 
realistic gravitational dual description. An appropriate set-up that suits the best for one's 
purpose should depend on the purpose. 

We will focus on qualitative aspects of hadron-virtual photon scattering amplitude at 
small X (at high center-of-mass energy). Since small x physics is dominated by gluon, not 
by quarks and anti-quarks, we do not find it a crucial element to implement flavor in the 
gravitational set up for the purpose of this article. For explicit calculation, we adopt the hard 
wall model [2], which is type IIB string theory onW x AdSs for some 5-dimensional manifold 
W with AdSs cut off at finite radius at infrared. Such a crude treatment of infrared geometry 
is sufficient for our qualitative study [2], and the choice of W becomes irrelevant (at least 
directly) for sufficiently small x [S]. Since it is almost straightforward to see how the AdS^ 
curvature and running of dilaton expectation value affects various observables in explicit 
calculations based on the hard wall model, one can also learn what happens in gravitational 
dual models that are asymptotically conformal or asymptotically free without carrying out 
calculations separately on these models. 

As an analogy of the electromagnetic global U(l) symmetry of QCD, we take a global 
symmetry of W in the gravitational dual. Since we are interested in the Compton tensoi§ of 
QCD, 

^(27r)V(p2 + g2-pi-gi)r'"^ = - [ I d\dSe-''^^-^e-^''^^-y{h{p2)\T{r{x)r{y)}\h{p{)), (1) 



we use the bulk-to-boundary propagator of an AdS^ vector field associated with a Killing 
vector of W in calculating the matrix element involving the global symmetry current. As 
for the target hadron in the gravity dual, we use a Kaluza-Klein state of a dilaton, whose 
wavefunction is given by a Bessel function in the hard wall model. Thus, the leading order 
contribution in l/N^ expansion is given by a closed string sphere amplitude with four NS-NS 



^In our convention, ry^i^ — diag( — , +, +, +). Let us remark that the Compton tensor T^'' in this letter is 
defined differently from one in [8]; the Lorentz indices ^, v are interchanged. 
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string vertex operator insertions |2]Jf| 

As we consider cases where the initial state "photon" or both the initial and final state 
"photons" are highly virtual, that is, qf ^ or g^,g2 ^ A^, the "photon" -hadron scatter- 
ing amplitude T^'^ can be decomposed into various contributions through operator product 
expansion of J^{x) and J'^{y) in QCD language. Such a decomposition still holds true in 
strongly coupled gauge theories (and hence in gravitational dual), except that the anomalous 
dimensions of operators in the expansion may be quite different from what one expects in 
the weak coupling regime. Reference [2] noted that the operators that are twist-2 in the 
weakly coupled regime still appear in the operator product expansion even in the strongly 
coupled regime, and their contributions to the Compton tensor T^'" dominate at sufficiently 
small X] this is because the "twist-2" contribution corresponds to exchange of leading Regge 
trajectory containing graviton in gravity dual language [9l [3]. We will thus focus on small x 
hadron- virtual photon scattering in gravity dual to study non-perturbative behavior of the 
"twist-2" contribution. 

Before writing down the Pomeron contribution to the scattering amplitude explicitly, let 
us note that the Compton tensor is described by five structure functions V^i,2,- -,5 as in [TU], 

T^-^ =V,P[q,rP[q,]; + V,{p- P^^ip ■ P[qi]r + ^3(^1 ■ P[q2]nq2 ■ PlqiW 

+ V,{q^ ■ P[q2]np ■ P[qi]r + V,{p ■ P[q2]nq2 ■ P[qi]r - Ae^^^f'^q^.q^,, (2) 

for a scalar target hadron, because of gauge invariance. In parity-preserving theory, A = 0. In 
the limit of purely forward scattering, the two structure functions of deep inelastic scattering 
are restored from Im Vi{x,r],t,q'^) — > g^) and (g^/(2x)) x Im V2{x,r],t,q'^) — )■ F2(x, g^). 

Here, we introduced a convenient notation 



P[q] 



(3) 



In this article, we will use the following notations. 



{Qi + (I2) 



{Pi + Q2] 



X 



2p ■ q 



9 • (gi - g2) 
2p ■ q 



(4) 



and t = — (gi — ^2)^ and s = = — (g + P^- 

^ The target hadron which is dual to a Kaluza-Klein state of a dilaton is a glueball. The case of a meson 
target can also be studied in the same way if we use open strings. For the case of a baryon target, wc should 
use ZJ-brane in the gravity dual. We will see that the saddle point value and singularities in the complex 
j-plane representation are important in describing the amplitude. Because they do not depend on the target 
hadron wavefunctions, they are expected to be unchanged even if the species of target hadron is replaced. 
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In the generalized Bjorken limit, A^, \t\ ^ qf, and for x much smaller than unity, the 
Pomeron contribution to the five structure functions are given by Jq and Ji |8] as in 

1 

V, - ^h, V, ^ ^h; (5) 

q^ 

Jo and Ji are given for vanishing skewedness rj in the form of 



For vanishing skewedness, the Pomeron kernel /C is p] 

1 f ,. l + e~'^^ 



sinvrj T^{j/2) 



as ^ 



the integration contour in the complex j-plane encircles the pole j = j^, and once the residue 
of this pole is picked up, a relation 

4 + z/2 

sets the (analytically continued) relation between spin j and anomalous dimension 7 = iv — j 
of "twist-2" operators in the large 't Hooft coupling A ^ 1 regime [3]. e^"^^^-* = {R/ z)"^ is the 
warp factor in the AdSs part of the metric in the hard wall model, 

rfs^UdSs = e2^(^)(r/^,da;''da;^ + {dzf), (9) 

and yj—g in ([6]) is that of this metric of 5-dimensional spacetime. R is the AdS radius, 
and the infrared cut off of the hard wall model at z = 1/A sets the confinement scale A. 
s = e~"^''^''e~'^^^'-'s, and a' is the slope parameter of the Type IIB string theory. '^[''J{t,z) in 
([7j) is the Pomeron wavefunction in the spin j channel, which is given by 



2i?sinh nu 



lUV^t/A) '''' ' \l i^UV^t/A) 

(10) 



^ More careful discussion on the choice of integration contour is given in [TTJ |S] . A pedagogical explanation 
of the origin of l/r'^{j/2) factor is also given in [8]. 
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in the hard wall modeljj 

The impact factor Phhiz') of the target hadron side is given by the normalizable mode 
wavefunction of the target hadron, as in Phh{z') = C0($(2;'))^. On the virtual "photon" side, 
the bulk-to-boundary propagator (non-normalizable wavefunction) of the graviton associated 
with the Killing vector of W is used; in the hard wall model, they are 

Pi,l\.{z) = c'jR'e-'^^^\{q,z){K,iq,zmq2z)K,iq,z)l (11) 

P?A^) = q2 Uz){K,{q,zmqlz)Ko{q2z)] (12) 

for Ji and Jq, respectively. k\ is a constant of a theory of mass dimension —3 and is pro- 
portional to A''^. 4, C0 and cj are dimensionless constants of order unity. See [5] for their 
definitions. 



3 Structure and Behavior of the AmpHtude 

3.1 Complex j-plane amplitude, Pomeron vertex and form factor 

Before discussing kinematical parameter (x, t,g^) dependence of the DDVCS amplitude in 
gravity dual, let us clarify a couple of conceptual issues associated with Pomerons. Using 
the explicit form of the Pomeron kernel ([7]) and Pomeron wavefunctions (fTOj) . amplitudes /j 
(i = 0, 1) in ([2]) can be rewritten (see [S] for details) as 



/i(a:, ?7 = 0,t,g^) ~ / dv 



l + e 



sinvrj^ J V'^{i^/2) 



[Ahh],. (13) 



where 



-1 I dz,rKz)P%{z)e-'^^^^ iW^^'^"^ 



1 



-tR) 



R. 

dz' -g{z')Phh{z') 



X (Rfi) 



(14) 



tz" 



(15) 



a parameter fi of mass dimension +1 is introduced in (IT^ [T^ in a way the observables 
are unaffected. One can change the integration variable of (fTS]l from u to j = j^; now the 
amplitudes /j are given by integration over the complex j-plane, and the contour becomes 
the one in Figure [1] (a). 



^Dirichlct boundary condition was imposed at the infrared boundary z = 1/A, just to make expressions 
simpler. 
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(a) (b) (c) (d) 

Figure 1: Singularities and integration contours in the complex j-plane. The hard-wall model 
is assumed for (a) negative t and (b) sufficiently large positive t, while a holographic model 
for asymptotic free running coupling is assumed in (c) with a smaller t, and in (d) with a 
larger t. Black dots are poles, wiggling lines in (a, b) are branch cuts, and open circles in 
(a-d) denote saddle points of the amplitude on the complex j-plane. 



The factor [C^'\j,q)] is now regarded as a function of j, and also depends on and 
but not on t or x. Its asymptotic form for ^ is given by 



1(3) 



2\j 



-C- 



(16) 



with a dimensionless constant of order unity cfj. that depends only on j. Here, 7(7) = iuj —j, 
and Uj = u{i) is the inverse function of j = jv x dependence and t dependence of the 
amplitudes /j come from the other factor [A/^/j]^. It can be rewritten as 



[A 



hh\ 



7(i) 



t/A) 



4\/A; 



X 



-) 



-t/A), 



(17) 



where 



t/A) is a dimensionless function of j and (v^/A). For the final expression. 



we used ~ q'^/x which holds at small x. Combining both, one finds that 



dj 



\di^j/d]\c 

r^(j72) 



-TTtJ 



sin(7rj) 



X 



A 



7(i) 



€^7ll,(v^/A). (18) 



This is in the form of inverse Mellin transformation, and the integration variable j is identified 
with the complex angular momentum (complex spin),§ 
Now, physical meaning of the separation between 



C) 



and [A, 



hh\ 



or 



<(v^))i 



IS 



clear. By changing the integration contour in the j-plane, (fT3| [TS!) can be rewritten as 



I ~ 



j€2N 



dj 



r'ij/2) 



1 



7(i) 



2V 



c)' 



-) 



7(i) 



t/A) 



(19) 



^Since we restrict ourselves to the scattering at 77 = 0, total derivative operators in field-tlieory language 
do not contribute to the OPE of the scattering amplitude. Thus, there is no subtleties in what this j is here. 
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This is regarded as an OPE form of Jj. The first factor in [■■■], which comes from ^C^^\j, q)]^, 
is regarded as the Wilson coefficient of OPE for a spin j G 2N operator; the parameter n 
is now identified with the renormahzation scale, because of its appropriate scaling behavior 
determined by the anomalous dimension 7(7) of the "twist-2" spin j operator. The second 
factor in [■ ■ ■ ] is identified with the spin j form factor, which is the coefficient of the [p'^^ ■ ■ ■ p'^^] 
term of the hadron matrix element of the spin j operator renormalized at the scale /j,. The 
gravity dual expression f[T^ justifies such an interpretation [T^ . 

Knowing physical meaning of these factors in the scattering amplitude f ll3p in a gravity 
dual model, one can define a GPD even in the model, which corresponds to a strongly coupled 
gauge theory. GPD as a function of x and t (we only consider the 77 = case in this letter) 
is defined as an inverse Mellin transform of form factors of twist-2 spin j operators (the 
second factor in [■ ■ ■] of (fT9|) ). The scattering amplitude /j is given by convolution of this 
GPD, inverse Mellin transform of the Wilson coefficient and that of the signature factor 
— [1 + e~'^*-']/ sin(7rj), just like in perturbative QCD factorization formula. The inverse Mellin 
transform of the signature factor gives rise to a light-cone singularity of a propagating parton 
(like the one in [13]), even in the gravity dual description. The GPD determined in this way 
is essentialljlll the same as Im /j, with of Im /j replaced by the renormahzation scale /x^; 
thus, various statements on Im /j in the rest of this section are also applied to the GPD after 
is replaced by /i^. 

Now that the field theory OPE interpretation of the gravity dual amplitude f|T3|l is clar- 
ified, let us go back to the amplitude (fT3|) and explicit expressions (iMl [T5|) once again. We 
will now clarify how this string theory amplitude on a warped background is related to the 
traditional Regge phenomenology ansatz. It should be noted that the DDVCS amplitudes 
Jj in gravity dual fll8p do not have a Pomeron pole like — a (t) ) in their j-plane repre- 
sentation apparently. There was once a pole l/(j — ju) at the stage of ([7]), but it is gone in 
( !T8|) . after picking the residue to evaluate an integral in ([7]). Nevertheless, one can see that 
the expression f lTSj) may have, in fact, many poles in the j-plane, rather than a single pole or 
none. 

To see this, we can use Kneser-Sommerfeld expansion of Bessel functions in the hard wall 

^Since GPD is defined as the inverse Mellin transform of form factors of "twist-2" spin j operators, it 
would become different when the normalization of the operators were changed in a j-dependent manner. 
We do not pay such a careful attention in this article. We claim similarity between Im Ii and GPD after 
replacement of by /x^ only at this level of precision. 
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model to rewrite [T'hhP*{j^t)]f. 



t/A) as [S| 



[ThhF*{j,t)]^ 
IhhFnU) 



E 

n=l 



■n(j) ^ — ) 



2fi 



7(i) 



'^5 



-g Phh{z) e 



-'^■^ ijl^\z) X {RAy. 



R' 



311/2 



(20) 
(21) 



The Pomeron trajectory (that contains graviton) of the Type IIB string theory on 10- 
dimensions (or on AdS^ after dimensional reduction on W) gives rise to a Kaluza-Klein 
tower of infinitely many Pomeron trajectories in hadron scattering on 3+1 dimensions. These 
trajectories are labeled by the Kaluza-Klein excitation level n; the masses of spin j G 2N 
hadrons are m, „, and their wavefunctions on AdS^ are tlji^\z). The factor l/{t — m?„) in 



( 120|) becomes a t dependent pole in the j-plane, the Pomeron pole, for any one of n's. In 
the hard wall model, the Pomeron trajectory (j = ap,n(^) relation set by t = r«|„) and the 
Pomeron wavefunction for the n-th trajectory are obtained holomorphically in j (not just for 
j G 2N) as in 



m 



R^ 



1/2 



(22) 



where, is the n-th zero of Bessel function 



rrij n is read out from the denominator in 



( ITSll : the wavefunction ipii [z] satisfies an equation of motion of a spin j field on ^^5*5, just 
like f irUj) does. Although explicit expressions above rely heavily on the hard wall model, con- 
ceptual understanding itself is quite general, and is applicable at least to any asymptotically 
conformal gravity dual models. 

Therefore, gravity dual descriptions of strongly coupled gauge theories come up with a 
following picture of Pomeron exchange amplitude. Individual Pomerons in the Kaluza-Klein 
tower couple to the target hadron with a coupling '~fhh¥n{j) in ^^M, which do not show any 
power-law fall-off behavior in large negative t. Only after all the Pomeron couplings '~fhh¥n{j) 
and Pomeron propagators l/(t — m|^„) oc l/(j — ap^„(t)) are combined as in (1201) . do we obtain 
what we might call a "Pomeron form factor" [T{j, t)]^, which has a power-law behavior in 
(see f l26|) ). Such a relation between a form factor of a conserved current and a combination of 
a Kaluza-Klein tower of hadrons, three point couplings and decay constants has been known 
for fixed spins (such as j = 1 and j = 2) [T3]. The relation fl20|) is regarded as an analytic 
continuation in j of the one for graviton (spin j = 2). 
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3.2 Saddle point in the j-plane 

It was a conventional wisdom of traditional Regge phenomenology that behavior of hadron 
scattering amplitudes at high energy are governed by the position of singularities in the 
complex j-plane. The same is true in gravity dual description of strongly coupled gauge 
theories. Singularities in scattering amplitudes in the complex j-plane representation depend 
on choice of gravity dual models. In case of the hadron- virtual "photon" scattering, however, 
the scattering amplitude can be approximated at saddle point in the j-plane (within a certain 
kinematical region which we call "saddle point phase" in §3.3p . In this case, the expression 
of the amplitude becomes not directly dependent on the singularities, and hence, detail of 
gravitational dual is irrelevent. In this subsection, we employ the hard wall model, and study 
the behavior of this scattering amplitude. 

In the hard wall model, there are no isolated poles in the complex j-plane for negative 
t — physical kinematics — except the branch cut that extends to negative j along the real axis. 
Figure [U^a); (t — m^^„) never vanishes for t < 0. Thus, the j integral of (fT8|) along the contour 
in Figure [T] (a) is evaluated by the saddle point method for small x [TT]. For < — t < 
the saddle point value of j is given bj 



v^ln(g/A) 



In 



(g/A)/(v^a:) 



(23) 



and 



Im /,(x, r/ = 0, t, g2) ~ ( - J di,, (v^/A) ; (24) 

the form factor g\^{^\l^-t j tC) has a dimensionless non-zero limit of order unity when —t — ?■ 0; 



it begins to fall off in power-lawo as 



. N -7(j)+2A-2 

^(v^/A) ^ ( ^ ) ~gl, (26) 



for larger momentum transfer A^ ^ —t. Here, A is the scaling dimension of the scalar field on 
AdS-^ containing the target hadron /i, and is a t-independent (but z/j-dependent) constant 



® The integrand of the Pomeron kernel ([7]) is rchable at |j| ~ 0(1), but not at \j\ > \f\ [9j. Therefore, 
we note that kinematical variables (x, and t) arc required to be consistent with j* ^ C(l)- 
^Thus, for ^ — t, the saddle point becomes 



^v* {q/K, X, ~t » A^) = VaM^^^. (25) 
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of order unity. Note, in particular, that the t-dependence of the scattering amphtude is given 
by the form factor that is once analytically continued to complex j-plane and then evaluated 
at the saddle point. The Regge factor (W^Y of string theory amplitude justifies focusing on a 
small range of j (or u) around the saddle point value at high-energy scattering; the power- law 
behavior in \/^ follows from the power-law wavefunction of the target hadron Ph^iz') and 
exponential cut-off of the Pomeron wavefunction, Ki^{^/^z) in f|T5l) in particular, in the 
limited range of Uj. 

The saddle point method provides a good approximation to the scattering amplitude for 
ln(l /x)/ -\/A ^ 1. It should be noted, however, that it allows us to keep all-order contributions 
in iu* = a/a ln(g/A)/ ln(l/x), which is not necessarily small and can be as large as 0{\^^'^). 
Thus, amplitudes and observables are expressed as functions of iu* (or j*). This makes easy 
to understand their dependence of kinematical variables (x,t, g^). 

Equation flMl) clearly shows the importance of the value of the saddle point of the j-plane 
amplitude. To see this more explicitly, let us define 

2n d\n[x Ii{x,r] = Q,t,q'^)] x / + 2n d\n[xli{x,ri = 0,t,q^)] 
^^^^"'^'^) = ' Ma:,t,g) = ^^^^ . (27) 

It is straightforward to see that they are given by 

7es(x, t, q^) = 7(j*), Aeff(x, t, g2) = f _ 1. (28) 

These effective exponents 7cff. and Aeg. depend on kinematical variables x, q^ and t only 
through the saddle point value j* . The ratio p = Re /j/Im /j = tan (|(j* — 1)) is also 
related directly to the saddle point value j* . 

We can see from (1251 [2^ that j* becomes large for large q^ and small for small x. Thus, 
at a given renormalization scale /i^ (replace q^ in Im /j), GPD in gravity dual still increases 
in the DGLAP evolution (that is, 7cfr < 0) for small enough x such that the saddle point 
value j* is still less than 2. Even at such a small value of a:, however, GPD eventually begins 
to decrease (that is, 7cs becomes positive) for large enough /i^. Such a behavior of GPD — 
qualitatively the same as in the real world QCD — in the DGLAP evolution was anticipated 
in [2]; this is indeed realized for finite yU.^ in gravity dual, when both q^ (/i^) and x dependence 
are included in the saddle point approximation. The other parameter Aefr characterizing the x 
evolution is known to increase gradually for larger q^ in the real world QCD [12] . As already 
seen in [TB], it does follow from gravity dual as well; we understand that this phenomenon 
is also essentially due to the increase of the saddle point value j* for larger q^. The same 
behavior is also obtained in perturbative QCD (See [TT]). 

The t dependence of the scattering amplitude is characterized by the slope parameter 
of the forward peak (also known as t-slope parameter), which we define for non-skewed 
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0.5 r 
0.4 
0.3 
0.2 
0.1 

12 3 4 

Figure 2: (color online) Slope B of the forward peak in DDVCS. The dimensionless value 
5 X is shown as a function of iu* f l23p : from top to bottom, blue (solid line) curve is 
for -y/— t/A ~ 0.01-0.1, red (long dashed) one for v^^/A = 1., yellow (dashed) one for 
\l^-tj = 3., and green (short dashed) one for ^/^/ = 6. For more details, see [8]. 

scattering as 

E,(x,r/ = 0,t,g2) = 2 |-lnImJ,(x,r/ = 0,t,g2). (29) 

at 

The t-dependence (and hence the slope parameter) comes entirely from the form factor for 
the physical kinematical region t < in the hard wall model. The t-slope parameter at t = 
in such a case can be regarded as the charge radius square of the hadron under "spin-j* 
probe" . Explicit expressions for the form factor in the hard wall model allow us to calculate 
the t-slope parameter; see Figure [21 The larger the spin j* (and hence ii^*), the smaller the 
slope. Therefore, through ( l23l) . the slope parameter decreases for larger g^, a prediction of a 
gravity dual model which cannot be made within perturbative QCD. 

We can also see from fl23| |25|) that the saddle point value j* depends weakly on ln(l/a;) 
or ln(iy^/A^) than on ln(g^/A^) for small x, and the dependence is in the opposite direction. 
Thus, the ln(l/a;) dependence (or ln(iy^/A^) dependence) of the slope parameter B must 
be weaker than its ln(g^/A^) dependence. This property of B, shared by 7eflf, Xcs and p, is 
an immediate consequence of the fact that the scattering amplitude is well approximated by 
the saddle point method on the j'-plane integral. This is a fairly robust feature of the saddle 
point approximation, and does not rely on specific details of the hard wall model. 

The saddle point turns out to be an important concept also in the scattering amplitude 
Im li in the impact parameter space, which is obtained by taking a Fourier transform in the 
transverse direction of the momentum transfer [pi — P2)- The x-dependent parton density 
profile in the transverse direction obtained in this way [18] in gravity dual shows Gaussian 
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profile at large impact parameter b, but is larger than the simple Gaussian form for smaller 
6 ([8]; see also [16]; deviation from the simple Gaussian profile is an immediate consequence 
of the fact that the 4D leading trajectory j = ap^n=i{t) is not perfectly linear). This core of 
larger parton density has approximately a linear exponential profile, e~'^J*'i''. The effective 
mass scale m^. „=i = rrin^i gradually changes as a function of b, and the linear exponential 
form smoothly turns into the Gaussian form for larger b, when iu* becomes of order unity. 
See [H] for more. 

3.3 Pole— Saddle Point Crossover 

Although we saw that the saddle point method well approximated the DDVCS amplitude 
for physical kinematical region t < in the hard-wall model, it does not in general. Even in 
small X, whether or not the scattering amplitude is well approximated by the saddle point 
method, depends on singularities of the amplitude in the j-plane representation, and hence 
on the gravity dual model one considers, and also the values of the kinematical variables x, 
and t. Although all the asymptotically conformal gravity dual models have a branch cut that 
is stretched to large negative j, there may also be some isolated poles in the j-plane as in 
Figured] (b). The hard wall model does not have such a pole for physical kinematical region 
t < (there are for sufficiently positive t), but there may be some for other UV conformal 
models that have different (and faithful to string theory construction) infrared geometry. 
Even more interesting are gravity dual models that are asymptotically free, where the cut is 
replaced by isolated singularities (Figure [T] (c, d)) [5]. 

When the saddle point (open circle in the figure) has a larger real part than any one of 
the singularities in the complex j-plane, then the integration contour in the j-plane should 
simply be chosen so that it passes through the saddle point, as in Figured] (a, c). When some 
of the singularities have larger real parts than the saddle point value j*, however, it is more 
convenient to take the contour as in Figure d](b, d), so that the scattering amplitude is given 
by contributions from finite number of isolated Pomeron poles j = an{t) (n = 1, 2, ■ ■ ■ ) and 
by a continuous integration over a contour passing through the saddle point. We refer to the 
two situations as saddle point phase and leading pole phase (or leading singularity phase), 
respectively. Such observables as Aes, 7eff, p and B exhibit totally different dependence on 
the kinematical variables x,q'^ and t in the two phases. In a given theory (i.e., in a given 
gravity dual model), one always enters into the saddle point phase for sufficiently large q or 
sufficiently negative t. In asymptotically free theories, it is likely that the leading singularity 
phase also exists for sufficiently small x and not so large negative t, even in the physical 
kinematical region t < 0. 
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The transition between the two phases is not singular but is a (smooth) crossover for 
finite X. This is because the saddle point approximation is never exact, and the "saddle 
point" should be thought of as a sort of diffuse object for finite x. Subleading singularities 
may also give rise to significant corrections to the amplitude simply given by the leading pole 
j = ai{t) for finite x, too. The transition becomes a singular phase transition only in the 
extreme small x limit. 



4 Lessons to Learn 



It is true that gravity dual calculation employs a background that corresponds to large 't 
Hooft coupling even at energy scale much larger than the hadronic scale A. Still, there are 
surprisingly many qualitative features in the gravity dual hadron-virtual "photon" scattering 
amplitude that are in common with the scattering amplitude in the real world QCD. Scatter- 
ing amplitudes in gravity dual have ln(g/A) and ln(l/x) scaling governed by 7cff = 7(j*) and 
-^cff = J* — 1 in the saddle point phase, and this is the same qualitatively as the prediction 
of the saddle point method in perturbative QCD, as we have already seen in §3.2[ The only 
difference between gravity dual and real world QCD is in the choice of anomalous dimension, 
7(j). Qualitative features are shared by both, and are controlled by the saddle point value 
3*- 

Qualitative features in t-dependence also show agreements. The gravity dual amplitude 
continues to the power-law fall-off behavior at large momentum transfer <^ —t. This 
property, which is expected to hold in the real world QCD theoretically [12] and confirmed 
experimentally, was difficult to be consistent with the traditional Regge phenomenology, but 
this problem is now overcome in gravity dual on warped spacetime (cf. [H [2])- Moreover, 
the t-slope parameter of (129!) and its result in Figure [2] for 77 = in gravity dual at saddle 
point phase nicely agrees with that in DVCS differential cross section |20], in that the slope 
parameter B decreases for larger ln(g/A), and is less sensitive to ln(l/x) or ln{W/A). Such 
observation suggests the (analytically continued) spin j form factors [Thh¥*{j,t)]fj, in both a 
gravity dual model and the real QCD are similar to each other. 

With so many basic qualitative features that gravity dual shares with the real world 
QCD, it is thus tempting to try to extract some lessons from the hadron-virtual "photon" 
amplitude in gravity dual. The origin of such similarity at the qualitative level becomes clear 
in the complex j-plane representation, where GPD is given by inverse Mellin transformation: 



H{x,v = 0,t;^^')r^ / A(± 



7(i) 



9iui 



t/A) 



(30) 
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Indeed, it is always possible to describe scattering amplitude by the j-plane integral in any 
theories, independent of whether the scattering is based on the real QCD or on the strongly 
coupled gauge theory studied in gravity dual; this is because Mellin transformation is only a 
mathematical transformation. This j'-plane integral also comes form OPE, notion of which is 
well-defined even in strongly coupled theories [2] . GPD in the j-plane representation (130|) is 
given by dropping the Wilson coefficient of OPE from the scattering amplitude Im/j f|T3|) . so 
the spin j form factor (reduced matrix element of twist-2 spin j operator) , which is the content 
of [■ ■ ■] in fl3Up . determines GPD. The spin j form factor is decomposed into two parts: RG 
evolution part {A/ fiy'^^\ and form factor at renormalization scale /i = A, g^j^^{y/^/ A); both 
show common properties in the real QCD and in gravity dual. The anomalous dimensions of 
the twist-2 spin j operator 7(j) in both theories are qualitatively similar [3], and g^^,{^J^-i/A) 
has the power-law fall-off behavior at —t ^ in common. 

The behavior of GPD is determined by the saddle point, or alternatively, by the leading 
singularity depending on which phase a set of parameters (x, t, yU^) sits in. This classification 
is applicable in any theories, not just in gravity dual. Then it is important to know which 
phase a given set {x,t,fi^) sits in. As we have pointed out, the behaviors of •jcs and Ace 
observed in HERA for DIS [15] are successfully explained by the predictions of the saddle 
point phase and are inconsistent with the prediction of the leading pole phase (or the leading 
singularity phase) [8]. Therefore, it is very likely that the (most of) kinematical region of 
DVCS that has been explored in HERA measurements is in the saddle point phaseo and 
GPD is approximately given by 

H{x,v = OAf^')- -] - gtJV^t/A). (31) 



X J \fi 



Most of the observed properties of the t-slope parameter B of DVCS in HERA [20] can be 
understood only from the fact that the kinematical region is in the saddle point phase (see 
[8]). A GPD model with a specific choice of gll ^^{\^^/ A) in [22] belongs to this category!"] 
One can also see that the saddle point expression f pT]) automatically satisfies a requirement 



10 In the standard parametrization of DVCS cross section '^'^dvcs(7-p^7p) _ ^ (l^)"' 

parameters {S, n) are given hy 5 = 4(j* — 1) and n = 7(j*) + 2j* in the saddle point phase. Thus, the saddle 
point phase implies rise of 5 for larger . HERA measurement [21] gives 5 = 0.44 ± 0.19 for = 2.4 GeV^, 
,5 = 0.52 ±0.09 for q^ = 3.2 GeV^, S = 0.75 ±0.17 for q^ = 6.2 GeV^, S = 0.84 ±0.18 for q^ = 9.9 GeV^, and 
,5 = 0.76 ±0.22 for q^ = 18 GeV^ in ZEUS, and S = 0.61 ± 0.10 ± 0.15 for q"^ = 8 GeV^, 6 = 0.61 ± 0.13 ± 0.15 
for q^ = 15.5 GeV2, and 5 = 0.90 ± 0.36 ± 0.27 for q^ = 25 GcV^, in HI. 

i^Reference [22] introduces an ansatz gf^^. ~ [j — a{t))~^{l — t/A'^)^^, inspired by a leading Pomeron pole 
(j — a{t))~^ and a power-law fall-off for (— t) A^. Our result ([20]) is conceptually different from this model; 
each Pomeron pole term with a Kaluza-Klein excitation level n does not show the behavior of power-law 
fall-off, but the power-law (|26p appears only after summing all the Kaluza-Klein tower of Pomeron pole terms. 



14 



that GPD should be consistent with DGLAP evolution, because /x-evolution is correctly taken 
into account in the j-plane expression ( l30l) . This is a nontrivial requirement on GPD modeling 
in general. One can consider, for example, a GPD profile given by PDF (GPD at t = 0) 
multiplied by some form factor at a given renormalization scale |23] : 




where B{x) = a'{l — x)'^ln(l/x) + ■ ■ ■ , and a' and p are parameters. The profile of GPD like 
this are not stable under DGLAP evolution. On the other hand, the GPD under the saddle 
point approximation ( 13T|) is given by the PDF multiplied by a spin j form factor evaluated 
at the saddle point value j = j*. The saddle point value j* depends on x and the factor- 
ization/renormalization scale fi. This result obviously takes into account renormalization 
effects, and hence is stable/reliable at any renormalization scale. 

The remaining task is to determine the spin j form factor at renormalized point fi = A, 
9iuj ci holomorphic function of j. This is along the line of the coUinear factorization 

approach (dual parametrization) to the modeling of GPD|2^. Derivation of g'jl^^ (-^/^/A) from 
the first principle is an impossible task in perturbative QCD, because of the non-perturbative 
origin of the form factor, and this is also hard in lattice simulation, because there is practically 
no way of finding analytic continuation of integer spin matrix elements into complex j. An 
alternative is to use predictions from the gauge/string duality, and a crude way is to use the 
prediction of the hard wall model derived in §3.1 as it is. Indeed, as we saw, the hard wall 
model can explain decreasing slope parameter B of DVCS for large q^, observed in HERA 
[20] . It is also possible to use more realistic gravity dual models for similar calculation, where 
at least we might want to require the model to have asymptotic free running for certain energy 
range (as in |25]) still with large 't Hooft coupling. 

If one wants to consider a gravity dual model that is truly dual to the real world QCD 
(if there is any), then it should run into a problem in its UV region of the geometry because 
of large curvature. This problem of gravitational description, however, may be alleviated 
by borrowing the understanding of perturbative QCD. Such strategy may not be totally 
nonsense. We saw that the singularities of the form factor in the j-plane are important in 
determining GPD, and gravity dual with asymptotic free running suggests that the singular- 
ities are infinitely many poleq^j [3]. The BFKL theory in perturbative QCD with a running 
coupling effect also suggests infinitely many poles in the j-plane [26]. Now, let us examine 

These poles correspond to trajectories of Kaluza-Klein modes in radial direction (z) of a single graviton 
trajectory in 10 dimensions (or on AdS^). On top of this tower structure, there is yet another tower structure 
of trajectories associated with the daughter trajectories of stringy excitations on 10 dimensions. 





X 




1 



(32) 
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how sensitive the position of the poles predicted from gravity dual are to the unreliable large 
curvature geometry in the UV region. In gravitational descriptions, each Pomeron pole has 
its wavefunction on the holographic coordinate, and the Pomeron wavefunction becomes lo- 
calized more and more into the IR region of the holographic radius when Re j of the pole 
increases. Therefore, the poles in large Re j are determined mainly by IR physics, and posi- 
tion of poles predicted by gravity dual should be reliable, while the poles in small Re j are 
quite sensitive to the unreliable geometry in the UV region. As for such smaller Re j poles, 
however, the position of the poles predicted by the BFKL theory (with asymptotically free 
running) will be reliable. Thus, by using both predictions from the gravity dual and the 
BFKL theory, the poles in the j-plane may be properly determined. 

In order to determine GPD completely, not only the position of the poles but also complete 
profile of the spin j form factor are required. The spin j form factor is given by integrating 
Pomeron wavefunction and impact factor in gravity dual, and in fact, also in the BFKL 
theory; the integration is carried out over the holographic radius z in gravity dual, whereas it 
is done over gluon transverse momentum kj_ in the BFKL theory. The similar structure in the 
k± factorization formula and the gravity dual scattering amplitude has been pointed out, 
and identification of gluon transverse momentum k± in the BFKL theory with holographic 
radius in the gravity dual is suggested [271 13]. Thus, one can retain the integration over 
the holographic radius in gravity dual in the IR (large z) region. The integration in the UV 
region may be replaced by that over k± coordinate in the BFKL theory; this large k± region 
is where perturbative QCD is reliable. 



Acknowledgments 

Part of this work was carried out during long term programs "Branes, Strings and Black 
Holes" at YITP, 2009 (TW), "Strings at the LHC and in the Early Universe" at KITP, 2010 
(TW), "High Energy Strong Interactions 2010" at YITP, 2010 (RN, TW) and also during a 
stay at Caltech theory group of TW. This work is supported by JSPS Research Fellowships 
for Young Scientists (RN), by WPI Initiative, MEXT, Japan (RN, TW) and National Science 
Foundation under Grant No. PHY05-51164 (TW). 

References 

[1] J. Polchinski and M. J. Strassler, Phys. Rev. Lett. 88 (2002) 031601, thep-th/0109l74| . 
[2] J. Polchinski and M. J. Strassler, JHEP 0305 (2003) 012, |hep-th/0209211| . 



16 



[3] R. C. Brower, J. Polchinski, M. J. Strassler, and C.-I. Tan, JHEP 0712 (2007) 005, 
[hi^h/0603115j . 

[4] HERMES Collaboration, A. Airapetian et al, Phys. Rev. Lett. 87 (2001) 182001, 
|liep-ex /0106068|; HERMES Collaboration, A. Airapetian et al, Phys.Rev. D75 
(2007) 011103, [liep-ex /0605108|; CLAS Collaboration, S. Stepanyan et a/., 
Phys.Rev.Lett. 87 (2001) 182002, f hi^x/0107043j ; CLAS Collaboration, S. Chen et 
al, Phys.Rev.Lett. 97 (2006) 072002, |hep-ex/06050l2] ; Jefferson Lab Hall A 
Collaboration and Hall A DVCS Collaboration, C. Munoz Camacho et al., 
Phys.Rev.Lett. 97 (2006) 262002, |n ucl-ex/0607029] ; Jefferson Lab Hall A 
Collaboration, M. Mazouz et al, Phys.Rev.Lett. 99 (2007) 242501, [0709.0450 
[nucl-ex]]; ZEUS Collaboration, S. Chekanov et al, Phys. Lett. B573 (2003) 46-62, 
P^ep^x/0305028] ; ZEUS Collaboration, S. Chekanov et al, JHEP 0905, 108 (2009) 
|arXiv:0812.2"5T7l [hep-ex]]: HI Collaboration, C. Adloff et al, Phys. Lett. B517 (2001) 
47-58, |hep-ex/0107005] ; and the three papers by HI in [20] • See also review articles [7]. 

[5] S. J. Brodsky, L. Frankfurt, J. F. Gunion, A. H. Mueller, M. Strikman, Phys. Rev. 
D50, 3134-3144 (1994) |hep-ph/9402"283] ; J. C. CoUins, L. Frankfurt, M. Strikman, 
Phys. Rev. D56, 2982-3006 (1997) |hep-ph/9611433 |. 

[6] D. Mueller, D. Robaschik, B. Geyer, F. M. Dittes, J. Horejsi, Fortsch. Phys. 42, 101 
(1994) |hep-ph/9812448]; X. -D. Ji, Phys. Rev. Lett. 78, 610-613 (1997) 



|hep-ph/9603249] ; X. -D. Ji, Phys. Rev. D55, 7114-7125 (1997) |hep-ph/9609381] ; 
A. V. Radyushkin, Phys. Rev. D56, 5524-5557 (1997) |hep-ph/9704207] . 

[7] X. -D. Ji, J. Phys. G G24, 1181-1205 (1998) |hep-ph/9807358] ; A. V. Radyushkin, 
|hep-ph/0101225| ; K. Goeke, M. V. Polyakov, and M. Vanderhaeghen, Prog. Part. 
Nucl. Phys. 47 (2001) 401-515, [hep-ph /0106012] ; M. Diehl, Phys.Rept. 388 (2003) 
41-277, |hep-ph/0307382 |; A. V. Belitsky and A. V. Radyushkin, Phys. Rept. 418 
(2005) 1-387, |hep-ph/0504030] ; S. Bofii and B. Pasquini, Riv. Nuovo Cim. 30 (2007) 
387, [0711.2625 [hep-ph]]. 

[8] R. Nishio and T. Watari, "High-energy Photon-Hadron Scattering in Holographic 
QCD," UT-11-02/IPMU11-0003/NSF-KITP-11-028, [1105.2999 [hep-ph]]. 

[9] S. S. Gubser, I. R. Klebanov, A. M. Polyakov, Nucl. Phys. B636, 99-114 (2002) 
[h ep-th/0204051| . 

[10] C. Marquet, C. Roiesnel, and S. Wallon, JHEP QA (2010) 051, [1002.0566 [hep-ph]]. 
[11] Y. Hatta, E. lancu, A. H. Mueller, JHEP 0801, 026 (2008) [0710.2148 [hep-th]]. 



17 



[12] E. Witten, Adv. Theor. Math. Phys. 2, 253-291 (1998) |hep-th/9802150] , S. S. Gubser, 
I. R. Klebanov, A. M. Polyakov, Phys. Lett. B428, 105-114 (1998) |hep-th/9802T09l , 
W. Mueck, K. S. Viswanathan, Phys. Rev. D58, 041901 (1998) |hep-th/9804035l . 

[13] The 3rd reference of [6]. 

[14] S. Hong, S. Yoon, and M. J. Strassler, JHEP 04 (2004) 046, [hep-th/ 03 12071]; 



S. Hong, S. Yoon, and M. J. Strassler, JHEP04 (2006) 003, |hep-th/0409ll8] ; S. Hong, 



S. Yoon, and M. J. Strassler, |hep-ph/0501197] . 



[15] E665 Collaboration, M. R. Adams et al, Phys. Rev. D54(1996) 3006-3056; ZEUS 
Collaboration, J. Breitweg et al, Eur. Phys. J. C7(1999) 609-630, |hep-e x/9809005]. 

[16] R. C. Brower, M. Djuric, I. Sarcevic, and C.-I. Tan, [1007.2259 [hep-ph]]. 

[17] R. D. Ball, S. Forte, Phys. Lett. B335,(1994) 77-86 [hep-ph/9405320] ; R. D. Ball, 
S. Forte, Phys. Lett. B336, (1994) 77-79. |hep-ph/9406385] . 

[18] M. Burkardt, Phys. Rev. D62, 071503 (2000) | hep-ph/0005 "l08], Phys. Rev. D 66, 
119903(E) (2002); J. P. Ralston and B. Pire, Phys. Rev. D 66, 111501 (2002) 
|hep-ph/0110075] ; M. Burkardt, Int. J. Mod. Phys. A18 (2003) 173-208, 
hep-ph/0207047 ]; M. Diehl, Eur. Phys. J. C25, 223-232 (2002), |hep-ph/0205208] . 



Eur. Phys. J. C31, 277(E) (2003). 
[19] G. P. Lepage, S. J. Brodsky, Phys. Rev. D22, 2157 (1980). 

[20] HI Collaboration, A. Aktas et al, Eur. Phys. J. C44 (2005) 1-11, fhep-ex/0506061]; 
HI Collaboration, F. D. Aaron et al, Phys. Lett. B659 (2008) 796-806, [0709.4114 
[hep-ex]]; HI Collaboration, F. D. Aaron et al, Phys. Lett. B681 (2009) 391-399, 
[0907.5289 [hep-ex]]. See also A. Freund, M. McDermott, M. Strikman, Phys. Rev. 



D67, 036001 (2003) |hep-ph/0208160| , and the references therein. 
[21] The 8th reference of [1] by ZEUS, and the 3rd reference of [20] by HI. 
[22] D. Mueller, |hep-ph/06050l3] . 



[23] M. Diehl, Nucl. Phys. Proc. Suppl. 161, 49 (2006) |hep-ph/ 05 10221] . 

[24] A. V. Belitsky, B. Geyer, D. Mueller, and A. Schafer, Phys. Lett. B421 (1998) 
312-318, |hep-ph/9710427] ; M. V. Polyakov, Nucl. Phys. B555 (1999) 231, 
|hep-p h/9809'483]; M. V. Polyakov and A. G. Shuvaev, |hep-ph/0207l53] ; D. Mueller, 



A. Schafer, Nucl. Phys. B739, 1-59 (2006) | |hep-ph/0509204] , K. Kumericki, 

D. Mueller, K. Passek-Kumericki, Nucl. Phys. B794, 244-323 (2008) f hep-ph/0703T79] . 

[25] O. Aharony, A. Fayyazuddin, and J. M. Maldacena, JHEP 9807 (1998) 013, 
|hep-th/9806159] . 



18 



[26] J. R. Forshaw, D. A. Ross, Quantum chromodynamics and the pomeron, Cambridge 
Lect. Notes Phys. 9, 1-248 (1997). 

[27] R. C. Brower, C.-I. Tan, Nucl. Phys. B662, 393-405 (2003) |hep-th/0207l44] ; 

S. J. Brodsky, G. F. de Teramond, Phys. Lett. B582, 211-221 (2004) |hep-th/ 03 10227] . 



19 



